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1 Introduction 

We consider here the problem of factorization of continuous matrix functions of the real variable. This means 
the representation of a given invertible square matrix G € (C(K.))” ” in the following form 

G{x) = G-{x)A{x)G+{x), (1.1) 

where continuous invertible matrices G~{x), G^{x) possess an analytic continuation in the lower 11“ = {z = 
X -\- iy : luiz < 0} and upper 11+ = {z = x -\- iy : Imz > 0} half-planes, respectively, and 

AM = <iiag((^) .. e„£Z, (1,2) 

The representation dm is called right (continuous or standard) factorization and can be considered for any 
oriented curve T of a certain classes which divides the complex plane into two domain D~ and 11+ with a 

( + \ 

J , S 11+, or for (if 0 S 11+). A similar representation 
G{x) = G+(x)A(x)G“ (x) 

is called left (continuous or standard) factorization. If the right- (left-) factorization exists, then the integer 
numbers asi,...,®„, called partial indices, are determined uniquely up to the order. In particular, there 
exist constant transformations of factors, such that aei > ... > as^. The factors G“, G+ are not unique. 
Relations of pairs of factors are described, e.g. in [53]. The right- (left-) factorization is called canonical 
factorization if all partial indices are equal to 0, i.e. aei = ... = = 0. 

Eactorization of matrix functions was first studied in relation to the vector-matrix Riemann (or Riemann- 
Hilbert) boundary value problem (see |20]1. The later was formulated by Riemann in his work on construction 
of complex differential equations with algebraic coefficients having a prescribed monodromy group (see, e.g., 
[H])- By using the method of the Cauchy type integral, the vector-matrix Riemann boundary value problem 
was reduced in [33], m to a system of the Fredholm integral equations. A part of the theory of the 
factorization problem is based on the study of such systems (see also [H]), though this approach does not 
answer, in particular, the questions of when it is possible to get factorization, how to construct factors and 
how to determine partial indices. 

Among other sources of interest to the factorization problems one can point out the vector-valued Wiener- 
Hopf equations on a half-line (see m, m, EZ]) and their discrete analogies, namely the block Toeplitz 
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equations (see, e.g., m. ED). The developed technique found several applications in diffraction theory, 
fracture mechanics, geophysics, hnancial mathematics etc. (see a brief description given, e.g. in [46] and 
references therein). 

Theoretical background for the study of the matrix factorization and its numerous generalizations is 
presented in [5], [T7], [33], HD (see also [33] )• The theory of the factorization is more or less complete (see 
0), but the above mentioned constructive questions about existence, factors and partial indices (which are 
very important for practical applications) have been answered only in a number of special cases. Among them 
one can mention rational matrix functions (see, e.g., |19jb functional commutative matrix functions (those 
satisfying G{t)G{s) = G(s)G(t), Vt, s G T, see [12])) upper- (lower-) triangular matrices with factorizable 
diagonal elements (see ED, EHD) certain classes of meromorphic matrix functions (see [I], 0, [30])) special 
cases of 2 X 2 Daniele-Khrapkov matrix functions (with a small degree of deviator polynomial) (see [15] . 
[I6] . [TT] . [28] 1. special cases of 2 x 2 matrix functions with three rationally independent entries (see [D, 
m, m and references therein), special cases of n x n generalization of the Daniele-Khrapkov matrix 
functions (see, e.g., [10], [26], [42], [46]), special classes of matrices possessing certain symmetry property 
(see [48] and references therein). Some approximate and asymptotic methods for matrix factorization have 
been developed too (see EH) ED) ED). For more detailed description of constructive approaches to solve 
factorization problem we refer to the recent survey by the authors [44] and references therein. 

In [85] we propose a new asymptotic method of construction of factors for a special class of n x n 
nonrational matrix functions in the case canonical factorization. The essential property of the considered 
matrices is that they become close (after suitable transformation) to a unit matrix (a similar assumption is 
used in El) 0). The idea to use such representation in factorization is related to that in general operator 
theory and has been exploited since the seminal work by Gohberg and Krein [^ (see also ED) ED). By 
reduction of the factorization problem to a special matrix boundary value problem we succeed in ED to get 
at each step of approximation a solution retaining the main properties of the given matrix. It helps to find 
conditions under which obtained asymptotic series converge. 

In the present paper we generalize the approach of [3^ and apply a new more general procedure in 
the case stable partial indices. Stability of partial indices is of special interest for factorization theory (see 
[6], [ID, and recent paper [2^) as it justifies the use of approximate procedures for factorization. In the 
considered case the stability of partial indices allows us to reduce auxiliary relations for components of our 
asymptotic formula to independent scalar Riemann boundary value problems. The crucial point is that 
some of these problems have negative index equal to —I. Thus, their vanishing at infinity solutions exist if 
and only if the given matrix function satisfies certain solvability conditions. Instead, we consider all above 
scalar Riemann boundary value problems in the class of functions bounded at infinity. It leads to corrections 
of the algorithm of [35]. Finally we get both factors of factorization problem in form of asymptotic series. 
The conditions of their convergence are found too. 

To the best of authors’ knowledge, our class does not coincide with any of the above mentioned classes. 
This class contains matrix functions which appears in the study of certain problems in fracture mechanics 
related to perturbation of the crack propagation (El) ED) ED - [40]). Another motivation is the use of such 
matrix functions in the study the inverse scattering problem (see [D). 

The paper is organized as follows. In Sec. 2 we introduce necessary notation and formulate the problem. 
A constructive algorithm is presented in Sec. 3. The method is illustrated by an example given in Sec. 4. 
We conclude our study by showing the quality of the factorization approximation by restricting ourselves 
only to the first asymptotic term and discuss the role of the chosen small parameter. 


2 A class of matrices. Problem formulation 

Let us introduce the following class of invertible continuous n x n,n > 2, matrix-functions SKn depending 
on a real parameter G R, satisfying the following conditions: 

(1) G^p G (C(K))”^” belongs to SKn if it can be represented in the form 

G^ = R^FR-\ (2.1) 

where bounded locally Holder-continuous on R (in general non-rational) invertible matrix R^ is such that 
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( 2 ) 


Rq — I, ( 2 - 2 ) 

(3) matrix function F does not depend on parameter has Holder-continuous entries fki G i.e. for 

all fc, Z = 1,..., n, 


\fki{xi) - fkiix 2 )\ < C 


Xl 


X2 


(2.3) 


and there exist 


lim fki{x) = fki{oo) = fki{-oo) = fkii+oo), 

X^OC 


satisfies the following asymptotic estimate at infinity 

(4) 

F{x) I, |x| —>■ oo, 

(5) F admits a right factorization, i.e. 

F{x) = F~{x)A{x)F~^{x), 


(2.4) 


(2.5) 


where Holder-continuous on M matrix-functions F (x), F'*'(x) possess an analytic continuation in the lower 
H” and the upper n+ half-plane, respectively, and A(a;) is a diagonal matrix 


A(a;) = diag 


x — i 


X + i 


X — i 


X + i 


aei,..., a3n G Z. 


(6) Partial indices asi > ... > ae„ are stable, i.e. 

aei - ae„ < 1. 

The matrices of the following form constitute a simple subclass of class SK 2 ' 


G^[x) = 


P{x) 

q{x)e 


q{x)e 


P{x) 


( 2 . 6 ) 


(2.7) 


In this particular case 


/ e’ 2 * 0 \ / p{x) q{x) \ 

Kix) = , F{x) = . (2.8) 

V 0 e-^ y V q{x) p(x) J 

This function appears after Fourier transforms of the Wiener-Hopf equation, describing a problem in fracture 
mechanics. 

We note that the matrix-functions of this type do not belong to any known class of matrix-functions 
which admit explicit factorization (see, e.g. [44] and references therein). 


3 An algorithm 

3.1 General construction 

By assumption any matrix Gi^(x) G SKn can be written in the form 

G^{x) = F-(x) (F-(x))~" G^ (F+(x))~" F+(x) =: F-(x)Gi,^F+(x), (3.1) 

where F~(x),F'^(x) are components of the factorization of the corresponding matrix F (see (12.51) 1. and the 
matrix Gi^^(x) is represented in the form (see (12.21) . (I2.4p l 

Gy^x) = (F-(x))"'g^ (i"+(a;))”' = (F-(x))-^R^FR-^ (F+(x))-\ 
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In addition to (l)-(6) we assume. 

(7) There exist a small parameter e = s(ip) (more exactly its value will be described later), such that for all 
a; € R and any finite (p 

Gi,^ix) = A{x) + eN^p{x), (3.2) 

with A(a:) determined by (12.51) and matrix N^p{x) being bounded and locally Holder continuous on M. 

Note that by assumption each entry of the matrix N^^ix) has a limit as |a;| —> +oo, i.e. there exists the 
value Nip{oo) = 0, i.e. Nip{x) € 77°(]R)lll Note also that the commutativity of the involved matrices is not 
assumed. 

Let us look for the first step in asymptotic factorization of the matrix Gi^ifi(x) in the form 

Gi,^(x) = (/ + eN-Jx)(A+(x))-^) A(x) (I + e(A-(x))-^N+^(x)) , (3.3) 


where 

A(a;) = A'''(a;)A“(a;), 


with 


A"*" (a;) = diag ■ 


X — i 

,x + i 


X — i 


(3.4) 

(3.5) 


where 

ae^ =max{0,8ej}, ae“ = min{0,sej}, j = l,...,n. (3-6) 

In stable case, when condition (1^ holds, we have two possibilities, namely, cBj > 0,Vj = 1, ... ,n, or 
EBa < 0,Vj = 1,... ,n. 

1°. Let us consider first the special case of matrix functions in SKn for which the partial indices of 
factorization (12.51) are such that 88^ = 1, Vj = 1,... , fc, eej = 0, Vj = fc + 1, ..., n, 1 < k < n. Denote by Ag 
the corresponding diagonal matrix in (12.51) 


Ao(a:) 


diag 



X — i 


, 1 ,... 


k times 


> . 


(3.7) 


In this case 

A(J'(x) = Ao(a:), Aq{x)=I. (3.8) 

In this notation representation (13.31) is equivalent to 

Gi,^{x) = (A+ (a;) + eA^i“,^(a;)) (J + eA^i+,^(a;)) . (3.9) 

Comparing terms for the lower powers of e we get the following matrix boundary value problem for 
determination of factors 

Ag {x)N+^[x) + N:[^^{x) = N^{x), a: G R. (3.10) 


It is customary to denote 

Mo^^{x) = N^{x). 

^As we can see later, the elements of our asymptotic factorization formula do not necessarily retain this 
property. 
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Since the matrix coefficient of problem (13.101) is a diagonal matrix, then this problem is equivalent to 
independent scalar boundary value problems 


(Si) 

nti{x) + n-^iix) = 


(ffl) 

ntnix) + = 


(ifi) 

+^fel(2;) = 

'^0,fcl7 

(S) 

= 

'^0,kn 7 

’^fc+l;l 

{x) + n'^+iaix) 

'^0,/c+l;l7 


(x) + n-f:^^.Jx) = 

^0,A:+l;n7 


^0,nl 7 

ntu{x) + n-^{x) 

^0,nn 7 


where n+g, are entries of matrices Mo^^p, respectively. 

Let us introduce a collection of pairs of analytic functions 


(3.11) 


^o[mQ^pg]{z) := 


OO 

z — i f mo^pq{t)dt 

2Tri J {t — i){t — z)' 

— OO 


z g (1 < P, 9 < n). 


(3.12) 


One can directly check that the introduced functions r2^[mo,pq](z) are analytic in the respective half-plane, 
and satisfy the identity 

Oq [?7ro,pg](x) Oq [77ZQ^pg](x) — 7TlQ pq(x), X g IR, (3.13) 

and the following additional condition holds true: 


[? 7 Zo,pq] (i) — 0 . 


(3.14) 


Using integrals (13.121) we get representation of bounded solutions to boundary value problems (13.111) . 

Bounded solutions to the first kxn boundary value problems can be delivered by the formulas (see, e.g. 
[201 P- 120]) 

= ^^^^0 1^0,pq]{z), n-g{z) =-%[mo^pq]{z), (3.15) 

for l<p<fc, l<g<n, and their boundary values are represented as 


npq{x) = ■^-^^^.{mo,pq{x) + So[rno,pq]{: 


(3.16) 


npq{x) = ^(mo,pq{x) - ^olm-o.pgKa::)), 

where Sq is a corrected singular integral operator on the real line (see, e.g. [101 P- 51-52]jl 


*S'o[mo,p9](x) 


OO 

x — i f mo^pq{t)dt 

J {t — i){t — x) 

— CO 


(3.17) 


We use the same notation for the Cauchy type operator and the singular integral operator S'o applied 
0 


to matrices: [Mo,,p](z), S'o[Mo,,p](a:). 
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On the other hand, bounded solutions to the remaining boundary value problems in (13.111) can be 
delivered by the formulas 


/c + l<p<n, l<g<n. 


and their boundary values are consequently represented as 

npq{x) = ^(mo,pq{x) ± 5'o[TOo.p9](a:)) ■ 


(3.18) 


(3.19) 


It is known (see, e.g. [301 P- 51-52]) that operator So is a bounded operator on the space Hp(R) of 
bounded Holder continuous functions on the open real line K. The norm of this operator will be denotecJl 

Cfi ■= (3.20) 

Let us refine the factorization of the matrix G'i,<p(a:), i.e. look for a presentation of Gi^ip{x) in the form 

Gi^^ix) = {1 + eiVf;^(x)(A+ (a:))"^ -b (a;))”^) x (3.21) 

XAo(x) (J -b eN+^{x) + e^N+^{x)) , 

where N^,^{x) are those found at the previous step. 

As before, relation (13.211) is equivalent to 

Ao{x) + eNq,{x) = (A+ (x) -b eN^,^{x) + e^N^^^{x)) (/ -b eN+^{x) + e^N+^{x)) . (3.22) 

Comparing terms at we get the following matrix boundary value problem for determination of 

factors 

= -N^^^{x)N+^[x), X e R, (3.23) 

where the right hand-side is already known. The solution to this problem is given by the formulas similar 
to (I3.15|) and (I3.18p . in which mo^pq should be replaced by mi^pq. Essential difference in comparison to the 
previous step is that the functions mi^pq(x) is no longer vanishing at infinity, but remain to be bounded. 
Anyway, the functions Vt^[mi^pq]{z) are well defined, i.e. analytic in the respective domains, having Holder 
continuous boundary values, and satisfy (13.131) . (I3.14|) . 

One can proceed in the same manner. Thus on the r-th step we use the representation 

Gi,<p(x) = Ao(x) -b eNq,{x) = 

= (/ + eNrJx){AUx))-^ + ... + e’-iV-^(x)(A+(x))-i) x (3.24) 

xAo(x) (/ -b eN+^{x) -b ... -b e^N+^{x)) , 

where N^^^{x),..., N~_^^^{x), N^^^{x),..., N^_i^^{x) are found in the previous steps. This leads to the 
matrix boundary value problem 


A^{x)N+ (x) + N Jx) = Alr-i^q,, X e: 


(3.25) 


where 


Mr- 1 ,q, = - [Afl,q,ix)N+_j^^^{x) + A^a^^ (x) 2 ,^( 2 ^) + ■ • ■ + ' 

The solution to this problem is given by a formulas similar to (13.151) and (13.181) . 

Finally, the factorization of the matrix function Gi^,^(x) is given in the form 

Gi,<^(x) = GC (x)Ao(x)GtJx), 


(3.26) 


^For the standard singular integral operator on the real line, i.e. the Hilbert transform, an exact values 
of its norm in Holder spaces are known, see [1]. 
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where the respective minus and plus matrix functions are presented by their asymptotic series 

CX) CO 

Gl,{x)=I + Y,e^N-J^){A+{x))-\ GtJx)=I + Y,e^N+Jx). (3.27) 

r=l r—1 

Here N-^{x), N+^ (x) are bounded solutions to the matrix problems (13.2511 for any r £ N. From the 
properties of the solutions to (I3.25P it follows that N~^{x), iV+^(x) are rather finite than vanishing at 
infinity (|x| —>■ +c»). 

2°. In general case of stable indices we have eej = s + 1,Vj = 1,... ,k, eej = s,yj = fe + 1,...,n, 
1 < k < n, where s £ Z is any integer number. In this case 


A(a;) 


diag 




V 

k times 


s+1 


X — i 

X + i 


S 



hence 

where Ao(a:) = A()'(a;) is the diagonal matrix given by dSH). 

Let us consider first an auxiliary relation 

Ao(x) + e ( 2 ^))”^) ^a{x) (/ + e7V+^(x)) . 

Then we arrive at the following matrix boundary value problem; 


Ao(x)IV+ (x) + iVi (x)= MoA^) = 


X — i 


N^p (x). 


(3.28) 


(3.29) 


Then, using the same line of the reasoning as for the previous case, one receives the asymptotic formula 
for factorization of the matrix-function (13.21) in the form 

Gi,^(x) = GZj^)A(x)GtJxl (3.30) 

where the respective minus and plus factors are presented in forms of asymptotic series 

00 00 

Gr,,(x)=H-^e’-iV-^(x)(A+(x))-\ G+Jx) = I + s’'■ (3-31) 

r—1 r=l 


3.2 Convergence 

The following Theorem gives conditions when, in the case of the stable partial indices, the proposed asymp¬ 
totic factorization becomes an explicit one, i.e. gives convergence conditions for the asymptotic series in¬ 
volved. 

Theorem 3.1. Let G^p be a matrix which meets the conditions (1) to (1). Let the parameter e (defined in 
(7)) satisfies the inequality 

|e| < IM (3-32) 

with the constant A = A(ip) being equal to 

A=\\NA)U1 + GA, (3.33) 
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||-^v’(')IU being the norm of the matrix function N^{x) in the Holder space and being the norm 

of the singular integral operator Sq : (K) —^ (R). 

Then both series in {3.27^ converge for all x G M. 

<1 It suffices to prove Theorem only in the above considered case 1° (see Subsec. 13.11) . It follows (see, 
e.g., [20l p. 48]) that singular integral operator S'o is bounded in Holder spaces since the “standard” singular 
integral operator S is. Let be the norm of So in Holder space H^{R). Then we have the following series 
of estimates 

||iVi^,^(-)IU < ai||^v(')IU(l + C^f), where ai = 

and 

l|Mi,^(-)IU<(ai||iV^(-)IU(l + ^^M)f: 

i.e. 

Il^ 2 ^v(-)IU ^ + C^)^ where Os = 

Finally, for each r > 2 

\\N?,4-)\\u<l\\Mr-iA-W + C^.), 

and 

||-^t--1,v3(')IIm — (o^lQ^i--! + a2Q;r-2 + . . . + ar-lCXi) ||./V,^(-)||” (1 + ^ , 

i.e. 

where Ur = ^(aia^-i + ■ ■. + Or-iai). 

Few first coefficients Ur we can calculate explicitly, namely, oi = 1/2, 02 = 1/8, = 1/16. As for 

coefficients with large enough indices we can proof by induction that 


OLf 


1 

16(r — 3) ’ 


Mr > 12. 


Therefore 


I < kl' 


(||iV^(•)IU(l + C^)")^ Vr>12. 


16(r-3)(l + C;,) 

Since the sequence ■^ i 6 (r- 3 )\i+c ) < 1 is increasing for sufficiently large r and 


lim 

r—^oo 


16(r-3)(l + C^) 


= 1 , 


then (since ||(A+(x)) ^|| = 1) the convergence of the series (|3.27l) for all x S R follows from (13.3211 . > 

Remark 3.1. It follows from the standard properties of the Cauchy type integral and singular integral with 
Cauchy kernel that conditions of Theorem \3.1\ guarantee convergence of the series in the right-hand side of 
1^3.27^ in the half-planes H”, n+, respectively. 

Remark 3.2. In fact, the decay of the second term in the right-hand side of i3.^) at infinity follows from 
the properties of matrices of the considered class and the proposed construction. Meanwhile, the solutions of 
matrix boundary problems i3.251) do not retain such behavior at infinity^ 

Remark 3.3. If the number A = A[ip) in Theorem \3.1\ is small enough, i.e. 

A=\\N^i-)Ul + C^r <1, (3.34) 


then the results remains valid for e 


1 and the described procedure will work without any changes. 


^This matrix norm can be any sub-multiplicative norm of square matrix-functions in iL^(R), i.e. satisfying 
an inequality ||i?(-) • C'(-)||^ < ||H(-)||^||C'(-)||^ for any pair of square matrix-functions B{x),C{x) (e.g. any 
induced operator norm, see e.g. [24] 1. 

®In contrast to the case of canonical factorization in [5^ , where the solutions of the corresponding problems 
at each step of factorization vanish at infinity. 


















3.2.1 The behavior of boundary values of solution at infinity 

Let us recall (see [20l § 4]) few facts on the behavior of the Cauchy type operator Hq and singular operator 
So, presenting these properties in the matrix form. 

1) Let n X n matrix function M € i.e. is Holder continuous on R and vanishing at infinity. The 

the Cauchy type integral n^[M]{z) (see (13.121) 1 determines two analytic matrices having boundary functions 
^l^[M]{x) € il°(]R) and satisfies the Sokhotsky-Plemelj formula 

^Iq [M]{x) — SIq [M](x) = M(x), X G M, (3.35) 

and vanishes at z = i: 

H+[M](i)=0. (3.36) 


Since 


OJIMK.) = T I 

— OO 

then it follows from (j3.35|) that there exists 


M(t)dt 


1 


M(t)dt 


{t — z) 2Tri J {t — i) 


lim n^[M]{z) = -n+[M]{i) = 0. (3.37) 

z—^oo 

The corresponding singular integral operator So[M]{x) (see (13.171) satisfy the relation 

lim So[M]{x) = -2n+[M]{i) = 0. (3.38) 

x—^oo 

It follows from (j3.37l) and another form of the Sokhotsky-Plemelj formula 

2Hq [M](x) = ±M(x) -I- So[M]{x). 

2) We can consider (13.351) as the boundary value problem (jump problem). Then any other bounded 
analytic solution Cl^[M]{z) to this problem (different from fl^[M]{z)) has the form 


H±[M](z) = H±[M](z) + H±[M](z), 


(3.39) 


i.e. differs of our special solution D,q[M]{z) in a constant matrix. The above considerations show that there 
exist 

lim Qq[M]{z) = —Oq [M](i), 


lim So[M]{z) = —2H){'[M](i). 

2^00 

Remark 3.4. Any solution to the problem iS. 1 fA) has the following behavior at infinity 

lim (z) = [^o.<p](d, 

Z—¥00 

where is a chosen bounded solution to the jump problem 

fl'^{x)+ fl~{x) = Mo,,p{x), X G K. 

If, in particular, = fl^[Mo^,p], then vanishes at infinity. 

The matrix Mx^,p{x) = —Nfi^(x)N^^(x) has the following asymptotics at infinity 

r 1 2 

lim Ml ^(x) = H(J'[Mo,(p](d , 

x—^oo L 

and is vanishing at infinity either matrix H()’[Mo,,^](i) or its square is the zero-matrix. 


(3.40) 

(3.41) 

(3.42) 


(3.43) 
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Remark 3.5. The algorithm described in Subsec. \3.1\ is not unique (cf. \35}/ ) since solution to the problem 
S3.25\) at each step of approximation is defined up to the constant matrix 


N+^ = [A+] ' - a-i} , N-^ = - - a-i} , 

where h^_i[Mr-i^ip\ is one of possible solutions to the jump problem 

+ VL~{x) = Mj—i^^p{x), a; € K, 


(3.44) 


and Cr-i is a constant matrix. 

The first k rows of are prescribed by equalities 

Cr — l,pq — 1 l,pq] (0 J ^ Si P Si k.^ S fj. q ^ n., 

but the remaining n — k rows can be chosen arbitrary. 

The choice made in Secs. \3.R \3.‘,A corresponds to 

Cr-i=n+_^[Mr-iJii), Vr = 0,l,2... (3.45) 

Remark 3.6. According to Remark \3.5[ is it possible to minimize the norm of the matrices at each 

consequent step (j = 1^2,... N) adopting the arbitrary constants. As a results, one obtains a factorization 
with required boundary behavior and asymptotic properties. The norms of the one or both factors 
can be controlled. However, this would not always guarantee the improvement for the convergence as the 
norm of matrix Mj^ip{x) to be normalized on the next step may well increase as the result. This suggests 
various scenarios for the optimization. One can also choose the arbitrary constant to preserve some local 
property of the solution (for example, to set a particular value of some components of the factors G^{ip,x) at 
a points x = Xg (depending on the number of arbitrary constants). Moreover, limiting values of the matrices 
at infinity (x = oo) also can be controlled in this way. All these fact allow to construct specific factorization 
with properties suitable for a particular application. 

Remark 3.7. Interestingly, whatever factorization is constructed, the factors satisfy the following conditions 

dpqi-i) = ^pq, ^<P<n, l<q<k; Gf^{oo) ■ G^(g(oo) = I. (3.46) 

Here Gf^{x) = { 9 ^ 79 ( 2 ^)}: Pi 9 = 1, 2,..., n and the integer k is defined in jS. 7| ). 

The authors believe that impossibility to set the uniqueness condition in the form Gf^{oo) = Gi^{oo) = I 
(as it was in the algorithm for the canonical factorization 135}/ ) is related to the fact that the proposed algo¬ 
rithms always produces the conditions automatically and thus it would impose too much restrictions 

to the factor if additionally Gf ^{ 00 ) = I is imposed. If any other algorithm for factorization of that class 
of the matrix functions is proposed which is free from the condition {3.4.&) i then probably some kind of a 
uniqueness result would be available. And finally, if one restricts the procedure to N asymptotic term only, 
the latter conditions will be valid with an accuracy of the order 0{ip^). 

All the above statements of this Remark are valid in the general case too (see i3.31\) } with the replacement 
of the matrices Gfg(a;) by Gf^{x). 


4 Example of a 2 x 2 nonrational matrix 

4.1 Asymptotic algorithm 

Example 1. Let us consider the matrix function on the real line R 

x—3i —2i 
x-\-i x-\-i 

42 x-\-3i 

x-\-i x-\-i 
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F{x) : 


(4.1) 










F{x) satisfies the properties (l)-(7) (Sec. [2]), i.e. F{x) € SK 2 in particular, it admits the right factorization 

F{x) = F~ {x)A{x)F'^ {x), X € R, 

with 

F-{x) = 


1 -1 


2 1 


-1 2 

Partial indices of this factorization are equal 


i^+(x) = I I , A(x) = Ao(x) = 


^ 0 

X-\-l 


0 1 


= 1 , = 0 , 


and thus they are stable *1 — *2 = 1- 
Then with R^{x) as in (12.81) 


G^{x) := R^{x)F{x)R^ {x) = 
By construction (see (EU) 


/ x-3i -2ie'‘f’^ 
x-\-i x-\-i 

x-\-3i 

x-\-i x-\-i 


Gi,^(x) := [F-{x)]-^G^ix) [i^+(x)]- 


Here (x)] ^ = F^{x), [T'+(x)] ^ = F (x). Thus 

Gi,<p(x) = 


( x—92+426 12i—4ie 

x+i x-j-i 


V 


-6i+4ie “^^+226"^^ x+9i—4ie 

x-\-i x-\-i 


As in general case, the matrix function Gi^tp{x) can be represented in the form 

Gi^(p{x) = Aq(x) + eN^. 

We can represent the matrix sN^p in the form involving trigonometric functions: 

eNp{x) = 


x-\-i 


Si sin^3 sin —2 cos^ \ 

fc+i 


—42 sin ^3 sin +2 cos 'j 
x-\-i 


16isin^ ^ 
ai+2 


/ 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


It follows from (|4.6D that the matrix function eN^p is “small”, i.e. vanishing at infinity with properly chosen 

e. 

Following Remark 13.31 we do not separate the existing small parameter e from the matrices but simply 
work the the entire matrices as they appear, i.e. consider representation 

Gi,vj(x) = Ao(x) + Mo^p{x). 

Now we give another representation of which is calculated directly 


Mo^p{x) = 


It can be factorized in the form 


! —8i+4ie +4ie*'^°^ 12i—4ie 

x-\-i x-\-i 

-62+426-^^^+226^^^^ 82-42e~^‘^^-42e^‘^^ 

CC +2 ai + 2 


Mo,pix) = M+^{x) + M^^^{x), 


(4.7) 
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where 


M+^{x) = 2* 


e-v 6_4e-'^»-2e~ 


MqJx) = 2i 


V 



x+i 

x+i 

3^gio>x|2g-o> 

4 2e*‘^*-2e-‘'’ 


x+i 

x+i 

2(e- 


2(e-'t>=:e-v>) 


X+2 

x+i 

2(e- 


2(g-i<4x g-o=) 


x-\-i 


x-\-i 


Let us find solution to the problem 

or what is equivalent to 4 separate scalar problems 


(4.8) 


(4.9) 


(4.10) 


x-i „ + 

x+t 'Ll 

x — i „ + 


nil = 




™ 0 , 12 , 

™ 0 , 21 , 

™ 0 , 22 , 


(4.11) 


where n^j, n^;, mo,fc/, A:, I = 1,2, are entries of the matrices respectively. We look for (partial) 

solution to the boundary value problems (|4.11l) in the class of bounded analytic functions. 

The bounded solutions to the problems 114.111) can be written in the form 

l(^) = (’^d.nW - cii) , n^i(z) = mo.ii + cii, 

liz) = (wdi2(2) - C 12 ) , n^2iz) = + ci2, 


{z) = - C21, 

(z) = m+{z) - C22, 


l(z) = mo2i +C21, 
,_{z) = mo22 + C22, 


where Cki,k,l = 1,2, are appropriate constants which serve to make n^i{z) to be analytic in 11+. 

Note that since the matrix M^^{x) tends to zero as x —>■ 00 , by construction, then, for any choice of 


the constant matrix Co, the matrix from the right hand side of (j3.2,8|) . i.e. 




has the following behavior at infinity 


Mi^,^(oo) = lim Mi^^{x) = Cq. 


(4.12) 


(4.13) 


We suggested before (see Remark [375|l to chose constant matrix Cq = {co,pq}p,q=i ,2 in such 
a way that 

Co = M„y w. 

It gives (Cj“’ = Co) 


Cr = >I'(v>)(3 J). 


(4,14) 


As the result, the limiting value of the new matrix-function Mi^(x) (see (I3.43p ) to be fac¬ 
torized on the next step (14.1311 is computed: Note that 




(4.15) 
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Note that according to the Remark [33], the constants C 21 , C 22 can be taken arbitrary. Just 
for illustration, we consider here the following three other choices of the matrix Cq\ 




4 -6 
0 0 


cf = «'(»’) 


4 -6 
0 4 


a 


(3) 




(4.16) 


Then, the corresponding limiting values of the respective matrix functions Mi^p{x) are 






Mfl(oo) = 0. (4.17) 


We chose to consider those cases for the following reasons: the hrst one looks the simplest 
natural choice (all arbitrary constants to set zero value). The second one was selected as it 
preserves the trigonal matrix form (with the equal elements on the main diagonal) for the 
matrices N^^{x) and thus for all consecutive matrices in the procedure. Finally, the third 

/Q\ 

case guarantees that the matrix function preserves the same property at inhnity as 

the initial matrix function 

The respective hrst order approximation of the factorization problem (I4.10p has the form 


Af.yyi) = (A„+)-' (M„y(i) - elf), Ne_^.j(x) = m^jx) + c, 


ki) 

-'0 1 


(4.18) 


In fact, we dehne here four different variants of the asymptotic expansion N^^{x) = 
N^^.j{x). Any of these solutions gives the hrst order approximates for the factorization in 
the form 

Ao(x) + Mo^^{x) = H{^^{x)Ao{x)H^^^{x) + 0 ( 9 ?^), (p ^ 0, (4.19) 

where 

/7,-^(i) = / + iV-^(i)(A+(a:))-‘. HtJx) = I + N+^(x), (4,20) 

and the estimate 0 (v?^) is uniform with respect to the variable x. 

For the next step of approximation we use representation similar to (I3.22p and arrive at 
the following matrix boundary value problem 


^0 + ^2,^) = = -N^^^{x)Nl^{x). (4.21) 

Entries of the matrix Mi^^p are calculated directly. 

By using solution to problem fl4.21l) we get approximation at the second step. 

A„(i) + Ar^(i) = ^ 0, (4.22) 

H^Jx) = I + iV,:^(a:)(A+(x))-‘ + Ni-Jx)(A+(x))-\ (4.23) 

H*Jx) = I + N+Jx) + N+Jx). 

The procedure can be naturally continued to an arbitrary order of the small parameter (p. 
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4.2 Numerical simulations 

4.2.1 Estimate of the reminder: standard choice of the constant Co = 

First, we consider the case with the choice fl4.14p of the constants, Cq^\ which guarantees 
the solvability of the problem flS.lOp and was instrumental to prove the convergency result. 
Having the first asymptotic terms of the factorization procedure in form fl4.20p . one can 
proceed to obtain the next consequent approximations (I4.23p by factorization of the matrix 
functions Mi^^p[x). In this specific example, one can write the analytical representation of 
Ml from fl4.21l) which, however, become cumbersome enough. For this reason we present 
numerical results for the components of the matrix Mi <^(a;) which highlight the fact that 
they are of the order 0{ip^) and uniformly bounded with respect to the variable x E R. 
On the other hand, the results can be treated as a measure of the factorization error if one 
decides to stop the asymptotic procedure immediately after the hrst step and to consider 
fl4.20p as the sought-for approximate solution. 

In Fig. [1] we present results concerning the two diagonal components, mjj{x, (f), j = 1, 2, 
of the matrix function Mi^^{x). Three different values of the parameter ip = 1,0.1 and 0.01 
are considered. Since, we expect that 

Mi,^{x) = p ^0, (4.24) 

the normalized values p~‘^mjj are given instead. Note that the estimate Ml^^p[oo) = 0{p^) 
follows immediately from (I3.43p and the choice of the hrst row in the matrix C^\ 




Figure 1: Diagonal elements, rrijj = mjj{x,p), j = 1 ,2, of the matrix defined in (14.121) 

for various values of the parameter p. The elements are normalised to the value of the parameter 
(p^. The horisontal lines show the limiting values of the normalised components at infinity 2(e“‘^ — 

Although the difference between the norms, ||mjj||oo, j = 1,2, and the limiting values of 
the components, \nijj{oo)\, differ less with the decrease of the parameter p (see Fig. [1]), the 
limiting value of the matrix function, Mi^^{oo), does not fully represent its max norm. For 
example, ||m22|loo >||^ii||oo apart from the fact that 1^22(00)! = |mii(oo)|. 

The other two component, rriij, i+j = 3, of the matric function Mi^^{oo) are presented in 
Fig. [2j To highlight the proven asymptotic estimate, here we aslo As expected, they vanish 
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Figure 2: The other two elements, rriij = mij{x,ip), i + j = 3, of the matrix Mi^^{x) for different 
values of the parameter tf = 1] 0.1; 0.01. It is clear that the both entries vanish at infinity rriij 0 
as Ixl —>■ oo. 


at infinity as x oo. One can also observe that ||mi 2 ||oo >||?^ 2 i|loo- Since, in the considered 
example, the functions are smooth at inhnity (meromorhc function?) one can expect that 

mij{x,(p) = 0{l/x), X ^ oo. (4-25) 

To demonstrate this fact, we also present in Fig. |3]the normalised values of these components, 
xip~‘^ ■ mij{ip,x), i + j = 3. Indeed the resulting functions are bounded and do not decay at 
inhnity. 




Figure 3: Non diagonal elements, mij(x), , f + j = 3, of the matrix Mi^^(x) for various values 
of the parameter (p. The elements are normalised differently to demonstrate their asymptotics at 
infinity. 

Surprisingly, their max norms grow while the value of the parameter (p decreases and, from 
the hrst glance, this contradicts to (I4.24p and (I4.25p . In fact, it is all true as the components 
mi 2 {(p,x) and m 2 i{(p,x) behave like sin^(<pa;) and thus both estimates (I4.24p and (I4.25p 
are valid uniformly with respect to one or the other independent variable. 
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4.2.2 Estimate of the reminder for the second choice of constant Co = Cq^^ 

Now we analyse the effect of a different choice of the constant matrix Co on the first step 
of the asymptotic procedure. Namely, we take it according to fl4.16P i which provides the 
second simplest option when we set the arbitrary constants 021,022 equal to zero. In Fig. 0] 
and Fig. 0] present similar results as in the Fig. [T]and Fig. [2] above. 

On the other hand, due to the structure of the chosen constant matrix Cq^\ the behaviour 
of the components in the first and the second lines are the same and their values are close 
to each other. 




X X 

Figure 4: Diagonal elements mjj{x) of the reminder matrix to be factorized on the next 

step for various values of the parameter ip. The elements are normalised to the value of the 
parameter ip^. 




X X 

Figure 5: The other two elements mij{x) of the reminder matrix Mi^^{x) for various values of the 
parameter ip. The elements are normalised to the value of the parameter 

A a direct consequence of this choice, the components m 2 i(c), m 22 {x) of the matrix 
vanish at inhnity (x —)■ 00 ) and behave similarly to the component mi 2 (x) and 
772,21 (x) from the previous example (compare fl4.25p ). It is also obvious from Fig. |4]and Fig. [5] 
that the estimate (14.241) is valid. However, and this was also suggested by the limiting value 
of the matrix Mi^^p[oo) in (14.171) . the max norm of the matrix function Mi^^(x) increase 
approximately in ten times in comparison with the same value from the previous example. 
Interestingly, the norm of the factors Nf^{x) in both those cases are the same. 
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( 2 ) 

4.2.3 Estimate of the reminder for the third choice of constant Cq = Cq 


It is important sometimes to have, if only possible, the factors well structured at inhnity and, 
here, a triangular matrix with equal diagonal elements comes as a desirable choice. This may 
be arranged by setting the constant matrix Cq in the asymptotic algorithm by the relation 
fl4.16P o. which guarantees solvability of the Hilbert-Riemann problem and, simultaneously, 
offers the both factors and the reminder represented by that type of the 

matrices at inhnity (compare fl4.16p 9 and fl4.17b l). 

In Fig. [6] and Fig. [7] we present results for this particular asymptotic expansion similar 
to those in Fig. [1] and Fig. [2j Now, apart of the fact we did not effect the max norms of 
the hrst asymptotic terms, the max norms of the reminder Mi ,^(x) increased even 

further what is the price for the chosen matrix structure. Thus the convergence may become 
slow and the convergence radius for the small parameter ip decreases on an order at least. 
As one can expect the maximal component of the reminder is m 22 - 




Figure 6: Diagonal elements mjj{x) of the reminder matrix to be factorised on the next 

step for various values of the parameter cp. The elements are normalised to the value of the 
parameter pP'. 


/Q'\ 

4.2.4 Estimate of the reminder for the last choice of constant Cq = Cq 

The last case we would like to discuss refers to the choice of the constant matrix Cq given in 
fld.lbp ci which immediately leads to fld.lZp q). This, in this particular case, the components of 
the reminder, Ml^^p{x), vanish at infinity similarly as it was for the original matrix function 
As the results, the second step of the factorization is repetition of the hrst one 
in terms of its practical implementation (e.i. there is no need to deal with the singular 
integral with the density having equal hnite values at inhnity (of the class Fr^(M)). Clearly, 
such choice of the constant Cq looks the most advantageous one from all the four examples 
considered in this paper. However, it is unlikely that such case is available for an arbitrary 
factorization. An interesting observation is that the hrst and the last examples are very 
similar in terms of the choice of the components of the constant Cq (compare fl4.14p and 
dlTbPa) and provide the most close values of the reminder Mi^^(x) in terms of the max 
norm. From comparing those examples, one can conclude that there is no much diherence in 
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Figure 7: The other two elements mij{x) of the reminder matrix for various values of the 

parameter ip. The elements are normalised to the value of the parameter ip^. 

terms of the asymptotic procedure to follow one or the others. Even if some decrease in the 
convergence rate may follow from a choice it may sense to proceed in this way if the desired 
property is guaranteed at the end of the procedure. Apart of the fact we can not prove this, 
judging from the bahaviour of the reminder to be factorized on the next stage, the standard 
choice (14.141) of the constant Cq provide the faster converging rate. The comparison also 
suggests that the minimization of the reminder norm is probably the most efficient strategy 
for the factorization. However, this restricts user’s change to deliver the factorization with 
some specihc properties of the factors. 



X X 

Figure 8: Diagonal elements mjj{x) of the reminder matrix Mi^^{x) to be factorised on the next 
step for various values of the parameter cp. The elements are normalised to the value of the 
parameter ip^. 

5 Conclusions 

We proposed the algorithm for effective asymptotic factorization of a class of nonrational 
matrix functions appeared in applications. We prove the convergence of the procedure 
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Figure 9: The other two elements mij{x) of the reminder matrix for various values of the 

parameter ip. The elements are normalised to the value of the parameter tp'^. 

thus, when the value of the small parameter lies in the convergence range, this factorization 
is explicitly represented in terms of series. The fact that the original unreturned matrix 
{(f = 0) possesses the stable set of the partial indices is crucial for the method. Although 
only right-hand factorization is considered here, the left-hand side analog follows trivially. 
The procedure does not lead to an unique form of the factorization factors and does not 
allow hxing factors limiting values at inhnity by a given constant matrix. However, the 
limiting values of the factors for any factorization are represented by the nonsingular mutually 
invertible matrices (see (Id.dbh o). We show effectiveness of the process on a numerical example 
where even the first step provides a very good accuracy while the next steps of the procedure 
are easy repeatable. Finally, by choosing arbitrary constants in a desirable way, one can 
direct the factorization process to preserve specific properties of the factors or to speed up 
the series convergence. 
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